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Electromagnetic modes are instrumental in building quantum machines. In this experiment, 
we introduce a method to manipulate these modes by effectively controlling their phase space. 
Preventing access to a single energy level, corresponding to a number of photons N, confined the 
dynamics of the field to levels 0 to A/" — 1. Under a resonant drive, the level occupation was found 
to oscillate in time, similarly to an N-\eve\ system. Performing a direct Wigner tomography of 
the field revealed its nonclassical features, including a Schrodinger cat-like state at half period in 
the evolution. This fine control of the field in its phase space may enable applications in quantum 
information and metrology. 


The manipulation of a quantum system usually in¬ 
volves the control of its Hamiltonian in time. An alterna¬ 
tive route consists in effectively tailoring its Hilbert space 
dynamically. This can be done by restricting the system 
evolution to a subset of possible states. When even a sin¬ 
gle energy level is disabled, the system evolution is deeply 
modified and is ruled by the so-called quantum Zeno dy¬ 
namics (QZD) [IHS]. As the name suggests, the level 
blockade can be realized by repeatedly checking whether 
the level is occupied, owing to the inherent back action of 
quantum measurements miais]. Alternatively, as in the 
present experiment, QZD can be achieved by blocking 
the level using a strong, active coupling to an ancillary 
quantum system [SHS], without any measurement [7]. 
These ideas have recently been demonstrated for atoms, 
using either Rb Bose-Einstein condensates [8] or Rydberg 
atoms [ 9 ]. However, the dynamics of these systems is in¬ 
trinsically confined to a finite number of energy levels. 
Here, using a circuit quantum electrodynamics architec¬ 
ture, we implement QZD of light. With its large number 
of energy levels and ease of control, a single electromag¬ 
netic mode offers a wider and more controllable phase 
space than atoms and two-level systems. 

The principle of our experiment is shown in Fig. 
One cavity mode of frequency fc is coupled to a qubit 
of frequency fq. For large enough detuning, their evo¬ 
lution can be described by the dispersive Hamiltonian 
hfca^a5-hfq\e){e\ — hxci^CL\e){e\, where h is Planck’s con¬ 
stant, is the ladder operator and |e) is the excited state 
of the qubit. The last term describes the frequency shift 
of the cavity (qubit) by —x, which occurs when the qubit 
(cavity) is excited by one extra quantum of energy. Ow¬ 
ing to this shift, a tone at frequency fq — Nx addresses 
only the transition between states \N) 0\g) and \N) O |e) 
for level widths smaller than x cni; here \g) is the ground 
state of the qubit. These levels then hybridize and repel 
each other. Their splitting is given by the Rabi frequency 
Qr at which the qubit population would oscillate in the 
case where the cavity is in state \N) (Fig.[^. Any transi- 
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tion to level N is now forbidden when the cavity is driven 
at resonance. Schematically, level N has been moved out 
of the harmonic ladder (Fig. [^. Then, starting from the 
ground state, the electromagnetic mode is confined to lev¬ 
els 0 to A" — 1, whereas the qubit remains in its ground 
state. The field dynamics is dramatically changed, resem¬ 
bling that of an A-level system, and nonclassical states 
similar to Schrodinger cat states develop. 



Figure 1: Principle of the experiment. Combined energy 
level diagram of the cavity coupled to the ancillary qubit. 
Each state is labelled as \k)^\g/e) for /c photons in the cavity 
and the qubit either in the ground (g) or excited (e) states. 
The level A = 3 is dynamically blocked by driving the qubit 
transition at /g — Ax (red wave). This drive hybridizes the 
states I A) (g) \g) and |A) (g) |e) and shifts their energy. A 
coherent excitation at the cavity frequency fc (green wave) 
induces a periodic evolution of the field confined to the first 
A levels. At half period, a nonclassical state similar to a 
Schrodinger cat is produced. 

In the experiment, we use the fundamental mode of 
a three-dimensional (3D) microwave cavity made out of 
bulk aluminum, which resonates at fc = 7.804 GHz. 
This mode is off-resonantly coupled to a superconduct¬ 
ing qubit El with bare frequency fq = 5.622 GHz and 
dispersive frequency shift x = 4.63 MHz. The cavity exit 
rate Xc = (1-3 @s)“^ is dominated by the coupling rate 
to two transmission lines connected to the cavity, which 
are used for both driving and the readout of the system. 
The relaxation rate 71 = (11.5 and decoherence 

rate 72 = (8.9 /is)“^ of the ancillary qubit are an order 
of magnitude smaller. 

The experiment is performed by first turning on the 
blocking tone at fq — Ax- For the level blockade to be 
effective, we choose CIr = 6.24 MHz much larger than 
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the level frequency width, which is about 7c. Then the 
cavity is driven at frequency fd ^ fc for a time t varying 
up to a few /is. The drive power is fixed throughout the 
experiment and would lead to an amplitude displacement 
rate Cd of about 3 in the cavity, were there neither 
damping nor nonlinearities. At time t, both the blocking 
signal and the cavity drive are turned off, and the field 
state is measured. Two measurement schemes are used 
to characterize the cavity state. Both methods use as a 
probe the same qubit that is used to provide the level 
blockade. 

The first method consists in measuring the probability 
Pk for the field to host k photons. To do so, a selective 
TT pulse is applied to the qubit at frequency fq — kx so 
that it gets excited if k photons are in the cavity. Mea¬ 
suring the probability to find the qubit in the excited 
state hence gives Pk m- This direct measurement relies 
on the qubit being initially in the ground state, which 
is not always the case. To account for this, two spuri¬ 
ous effects need to be considered. First, the qubit has 
a residual thermal population of 22 %. Second, imper¬ 
fections of the blockade induce a parasitic excitation of 
the qubit, which was measured in time and remains be¬ 
low 18 %. Taking these effects into account leads to an 
accurate determination of the conditional probability 
of measuring k photons, given the qubit being initially in 
the ground state citeFootnote2. 

The resulting probabilities are shown in Fig. as a 
function of time for several photon numbers k and for N 
from 2 to 5. Levels with more than N photons are unoc¬ 
cupied, as expected from Hilbert space blockade. Early 
in the evolution, the occupation of the levels k > 1 rises 
in the order of increasing energy, similarly to a coherent 
state of increasing amplitude (Fig. S3c in [l3]). Later, 
the level distribution ’’bounces” off a wall at /c = A^, so 
that the probabilities start to oscillate. The period of 
the oscillations increases with N as expected, because it 
takes more time to reach N — 1 photons with a constant 
drive as N increases. The case = 2 is straightforward 
as it implements an effective qubit [8]. The time traces 
of Fig. correspond to Rabi oscillations of a two-level 
system. For larger the evolution is similar to that 
of a resonantly driven A/"-level system, as seen in Ryd¬ 
berg atoms [9]. In particular, at half period, Pq displays 
plateaus all the more pronounced as N gets larger. Fi¬ 
nally, the A^ — 1 level occupation evolves in opposition to 
level 0, with a maximum at half period. 

The evolution of the field can be modeled with the 
Hamiltonian of the cavity 

H = ihed{a\^ - aAr)/27r + h{fc - fd)aN'^N - hX{a\^)‘^a%, 

( 1 ) 

written in the frame rotating at the drive frequency fd- 
These three terms describe the coherent drive, its detun¬ 
ing, and the nonlinearity of the cavity, respectively. The 
blockade at level N suppresses transitions to and from 
\N). This is entirely modeled by an effective annihila- 
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Figure 2: Evolution of the photon number probabil¬ 
ities. Measured (dots) and theoretical (solid lines) photon 
number state probabilities Pk as a function of time t. The 
blocked level N ranges from 2 to 5 and is indicated on each 
panel. The standard deviation, estimated from the curve 
FW(t), is below 1.7%. 


tion operator 

ajv = a-\/]V|Af-l)(Ar|-VivTT|Ar)(Ar+l|, (2) 

for which these forbidden transitions are removed. Fi¬ 
nally, the relaxation at rate Xc is modeled using a master 
equation of Lindblad form m- 

Using this model, one can reproduce with good agree¬ 
ment the photon number probabilities Pk (Fig. [^, with 
the displacement rate the drive detuning fc — fd^ aiid 
the cavity anharmonicity A as the fit parameters. The 
displacement rate which directly determines the pe¬ 
riod of oscillations in the level occupation (Fig. [^, is 
found to be finely tuned between 2.83 and 3.05 /is^ for 
each blocked level N. The anharmonicity A = — 70 kHz 
is negative, as confirmed by independent measurements 
(Fig. S2 in [IS]). This negativity is in contradiction to 
all reported transmon qubit and 3D transmon models 
such as those in Refs. [HHisi. At last, we attribute the 
origin of the detuning fc — fd to an energy shift of the 
levels caused by the blocking tone. It was found to be 
—0.1 MHz for A^ > 2 and —0.4 MHz for N = 2. These 
values are consistent with the system being more strongly 
disturbed by the blocking field at A^ = 2 than for larger 
N. 

These measurements demonstrate how the cavity is 
transformed into an N-levei system by dynamically in¬ 
ducing a blockade at an arbitrary level N. A full charac¬ 
terization of the field, which goes beyond measuring only 
photon number probability, requires tomography. The 
Wigner function is a complete representation of the quan¬ 
tum state of the field in continuous variables. It can be 
expressed as W{a) = {DaPDj^), where ^ is 

the field displacement operator and V = the pho¬ 

ton parity operator. A direct Wigner tomography of the 
field is thus performed by first displacing it by a complex 
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Figure 3: Time evolution of the Wigner function of a quantum mode of light under QZD. Measured (top rows) 
and predicted (bottom rows) Wigner functions lF(a) as a function of the displacement amplitude o, for a blockade at iV = 3 
(top panel) and N = 4 (bottom panel). The time t (nonlinear scale) for the frames shown in A to I is given above each panel. 
The field is confined in phase space by a barrier at amplitude |o| = VN (white dashed circle). Negative values of the Wigner 
function, in blue, demonstrate the nonclassical nature of the field produced under QZD. The model used is the same as for 
Fig. Wigner functions are directly measured here and not reconstructed. 


amplitude a and then measuring the average photon par¬ 
ity. This can be realized by mapping even and odd pho¬ 
ton numbers on the qubit ground and excited states [13^ 
HH. In practice, this is performed in three steps. The 
qubit is first prepared in state (|^) + |e))/\/2 using a 
pulse at the qubit frequency. The 7r/2 pulse bandwidth 
is designed to be larger than 10% so as to succeed regard¬ 
less of the cavity state. Then, the qubit is left to evolve 
freely during a time r, where it acquires a phase shift 
that depends on the cavity state owing to the evolution 
operator e^^Trxra a|e)(e| waiting time is chosen to be 

r = 1/2% = 108 ns, so that the qubit phase increases 
by kir when the cavity state is \k) (Fig. S4 in [l3]). Fi¬ 
nally, another broadband 7r/2 pulse is applied in order to 
flip the qubit to the excited (ground) state in case of an 
even (odd) number of photons. Measuring the qubit state 
then directly leads to the Wigner function in contrast 
to indirect reconstruction procedures such as maximum 
likelihood. 

The measured Wigner functions for the field state are 
shown in Fig. at various times during the first oscil¬ 
lation period for blocked levels N = 3 and N = 4. At 
time zero (Fig. Ih), the cavity is in the vacuum state 
and the Wigner function is a positive Gaussian distri¬ 
bution centered at a = 0, with a width reflecting zero 


point fluctuations. As time increases, the resonant drive 
displaces the Wigner function along the real and posi¬ 
tive axis. Without the blockade, and neglecting losses 
and anharmonicity, the vacuum state observed at time 0 
would simply get displaced hy a = Cdt. In the presence 
of the blockade at level A, all levels with k > N remain 
empty so that the field cannot be in a coherent state, 
and the distribution gets distorted. Indeed in Fig. |^, 
the quasi-probability distribution seems to hit a wall in 
phase space at an amplitude There, it undergoes 

a rapid counterclockwise evolution (Fig. to H) along 
the corresponding barrier (white dashed circle). After a 
full oscillation period (Fig. |^), the cavity goes back to a 
state close to the vacuum state of Fig. Ik- These snap¬ 
shots of the field state are close to the ones predicted 
in Ref [5] for light under QZD, and similar to the one 
observed in the levels of a Rydberg atom [9]. 

Besides confirming the confined and periodic evolution 
of the field under QZD, this tomography reveals the for¬ 
mation of nonclassical field states, as indicated by the 
appearance of negative values in the Wigner function. 
These negative values develop while the field undergoes 
a transition along the barrier in phase space. At half 
period, the state is close to an equal superposition of 
positive and negative field amplitudes at 4z\/N — 1. The 
quantum nature of this superposition manifests itself by 
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Figure 4: Wigner tomography at half period. Measured 
(top row) and calculated (bottom row) Wigner functions of 
the cavity field for various blockade levels N from 2 to 5, 
taken at half period (t = 0.51 /xs, 0.64 /xs, 0.7 /xs and 0.75 /xs). 
The color scale is rescaled compared with Fig. by A = 0.7. 
Similarly to Schrodinger cat states, these states exhibit fringes 
with alternating positive and negative values. 

the appearance of interference fringes in the Wigner func¬ 
tion. These characteristic patterns can be compared for 
all blockade levels 2 < < 5 in Fig. The photon par¬ 

ity is that of the highest allowed level N — 1^ which is also 
the number of negative fringes. These properties are rem¬ 
iniscent of a quantum superposition between two coher¬ 
ent states of amplitude so-called Schrodinger 

cat state m- The evolution of the Wigner function can 
be predicted using the model above. The calculated dis¬ 


tributions are shown below each measurement frame in 
Figs. [^and|^ The observed asymmetry between positive 
and negative imaginary amplitudes is attributed to drive 
detuning, losses and anharmonicity nuEU]. The model 
reproduces qualitatively the measured Wigner functions 
for all times and values of N. 

In conclusion, we have demonstrated Quantum Zeno 
Dynamics of microwave light. Our experiment shows 
that a mode that cannot contain N photons undergoes a 
coherent evolution that is confined in phase space. Effec¬ 
tively, an electromagnetic mode is transformed into an 
N/2 spin. As a result of this transformation, exotic non- 
classical states develop under a simple resonant drive and 
share similar features with Schrodinger cat states of light. 
Our work could be extended to allow evolution outside 
of the exclusion circle, where the barrier acts as a pe¬ 
culiar scatterer in phase space [5]. Other cat-like states 
and squeezed states could then be produced. Tailoring 
the Hilbert space in time is an important new resource 
for quantum control. Our experiment enables the real¬ 
ization of various protocols based on QZD, such as gener¬ 
ation and protection of entanglement [2TII23] . and quan¬ 
tum logic operations [24]. Moreover, by turning on and 
off in time the blockade level and combining it with fast 
amplitude displacements, which is feasible by improving 
the cavity decay rate by an order of magnitude, it would 
be possible to realize phase space tweezers for light mu- 
A possible and important outcome would then be the ma¬ 
nipulation of Schrodinger cat states in a unique way and 
the quantum error correction of cat-qubits, which are a 
promising quantum computing paradigm |25| . 
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